February 1, 2008 9:50 WSPC/INSTRUCTION 

IWARA'NS'2005Tublished 



FILE 



International Journal of Modern Physics D 
© World Scientific Publishing Company 



A NOTE ON THE LARGE- ANGLE ANISOTROPIES 
IN THE WMAP CUT-SKY MAPS 



A. BERNUI 

Instituto Nacional de Pesquisas Espaciais - Divisdo de Astrofisica 
Av. dos Astronautas 1758 
12227-010 Sao Jose dos Campos - SP, Brazil 

B. MOTA, M.J. REBOUCAS 

Centro Brasileiro de Pesquisas Fisicas 

Rua Dr. Xavier Sigaud 150 
22290-180 Rio de Janeiro - RJ, Brazil 

R. TAVAKOL 

Astronomy Unit - School of Mathematical Sciences 
Queen Mary, University of London 
Mile End Road, London El 4-NS, UK 

Received Day Month Year 

Recent analyses of the WMAP data seem to indicate the possible presence of large- 
angle anisotropy in the Universe. If confirmed, these can have important consequences 
for our understanding of the Universe. A number of attempts have recently been made 
to establish the reality and nature of such anisotropics in the CMB data. Among these 
is a directional indicator recently proposed by the authors. A distinctive feature of this 
indicator is that it can be used to generate a sky map of the large-scale anisotropies of the 
CMB maps. Applying this indicator to full-sky temperature maps we found a statistically 
significant preferred direction. The full-sky maps used in these analyses are known to 
have residual foreground contamination as well as complicated noise properties. Thus, 
here we performed the same analysis for a map where regions with high foreground 
contamination were removed. We find that the main feature of the full-sky analysis, 
namely the presence of a significant axis of asymmetry, is robust with respect to this 
masking procedure. Other subtler anomalies of the full-sky are on the other hand no 
longer present. 

Keywords: Observational cosmology; cosmic microwave background; large-scale 
anisotropies in CMB; large-angle anomalies in CMB. 



1. Introduction 

The wealth of high resolution data provided by the Wilkinson Microwave Anisotropy 
Probe (WMAP has confirmed to very good approximation the standard cosmo- 
logical picture, which predicts a statistically isotropic Gaussian random cosmic mi- 
crowave background (CMB) temperature fluctuations. Despite this success, several 
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large-scale anomalies in the CMB have been reported including indications of non- 
Gaussianityp2~El evidences for a North-South asymmetrypl and the so-called "low-£ 
anomalies" such as the surprisingly small values of the CMB quadrupole and oc- 
topole moments^ and the alignment of the quadrupole and octupole momentsp"^ 
(in this connection see Ref . [T0|) whose direction has been suggested to extend to the 
higher multipolejm ( se e also Ref. [12] for a detailed discussion) . In addition, there 
are also in dication s for a preferred axis of symmetry or directions of maximum 
asymmetryff^Hin] 

The possible origins of such unexpected anomalous features of CMB are at 
present the object of intense investigation, with several potential ex plan ations, in- 
cluding unsubtrac ted foreground contamination and/or systematics unconsid- 
ered local effects, ^^2 other mechanisms to break statistical isotropyj^ and also 
extra-galactic origin (see Refs. 13 Qj] and IT3HI51 for details, and Ref. [21] for re- 
cent related references). If they turn out to have a cosmological nature, however, 
they could have far reaching consequences for our understanding of the Universe, 
in particular for the above-mentioned standard cosmological scenario. 

Recently we proposed^ a new dir ecti onal indicator a — cr(9, 4>), based on pair 
angular separation histogram (PASH)j22lto measure large- angle anisotropy in the 
WMAP data. An important feature of our indicator is that it can be used to gen- 
erate a sky map of large-angles anisotropies from CMB temperature fluctuations 
maps . W e have produced and studied in details ct— maps generated from the full-sky 
LILCPS 'cleaned' TOH, 8 and co-addecP WMAP maps, and found a statistically 
significant preferred direction in th ese W MAP maps, which agrees with the preferred 
asymmetry axes recently reported! ^^ These results were found to be robust with 
respect to the choice of the full-sky WMAP CMB maps employed. However, since 
full-sky maps are known to have residual foreground contamination^! and c qmpli- 
cated noise properties P their choice in the "low-£" studies is not a consensus OSES] 
Thus, the question arises as to whether our results hold for cut-sky maps. Our main 
aim here, which extends and complements our previous workp^l is to address this 
question by considering the LILC map with a Kp2 sky cut. To this end, in the next 
section we give an account of our large-angle anisotropy indicator, while in the last 
section we apply our indicator to the LILC map with a Kp2 sky cut, and present 
our main results and conclusions. 



2. Large-angle Anisotropy Indicator 

For a detailed discussion of the indicator briefly presented in this section we refer 
the readers to Ref. l22l 

The key point in the construction of our indicator is that a homogeneous distri- 
bution of points on a two-sphere S 2 is seen by an observer at the center of S 2 as 
isotropically distributed, and therefore deviations from homogeneity in this distri- 
bution give rise to anisotropies for this observer. 

Mutatis mutandis, since in CMB studies the celestial sphere is discretized into a 
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set of equal size pixels, with a temperature fluctuation associated to each pixel, the 
idea in the CMB context is then to construct an indicator that measures deviation 
from homogeneity in the distribution of pixels with similar temperature. The first 
step towards the construction of this indicator is subdivide a given CMB map 
into a number of submaps, each consisting of equal number of pixels with similar 
temperatures. The next step is to devise an indicator to measure the deviation 
from a homogeneous distribution of these pixels. The construction of our indicator, 
cr = (i(9, 4>) is based on angular separation histograms (PASH), which are obtained 
by counting the number of pairs of pixels whose angular separation a lies within 
small sub-intervals (bins) Ji € (0, 7r], of length 8a = ir/Nbms, where 



. 8a 8a 
■A = a>i , a, + — - 



i = l,2 ) ...,JV, 



bins 



with the bin centers at a, = (i — h) 8a . The PASH is then defined as the following 
normalized function 

= 7 ~ 7T j~ E V(<») , (2) 
n(n — 1 ) 8a z — ' 

where n is the total number of pixels in the submap, rj(a) is the number of 
pairs of pixels with separation a 6 J 2 :, and where the normalization condition 
* $(q^) Sa = 1 holds. 
Now, for a distribution of n pixels in the sky sphere S 2 one can compute the ex- 
pected number of pairs, rj exp (a) with angular separation a £ Jj. From this quantity 
one obtains the normalized expected pair angular separation histogram (EPASH), 
which is clearly given by 

$ex P {ai) = — — E Vex P (a) = P{ai) = V{oti) , (3) 

aeJi 

where N = n(n — l)/2 is the total number of pairs of pixels, P(ai) = 
X)aGJ Vexpicx) / N is the probability that a pair of objects can be separated by an 
angular distance that lies in the interval Ji, V(a,i) is the corresponding probability 
density, and where the coefficient of the summation is a normalization factor. Equa- 
tion ([3]) makes it clear that the EPASH Q e xp(a) gives the distribution of probability 
of finding pairs of points on the sky sphere with any angular separation ai £ (0, 7r]p 
We denote the difference between the mean PASH (MPASH), ( $ b s (ai) ), calcu- 
lated from the observational data, and the EPASH Q exp (ai), obtained from an 
statistically isotropic distribution of pixels, as 

T(oi) = ( $ o6s (a 4 ) ) - ^expioti) . (4) 



a For a homogeneous and continuous distribution of points on S 2 all angular separati ons < a < it 
are allowed, and the corresponding probability distribution can be calculated to givd^SJ & exp (ct) = 
■p (a) = i sin a . This is the limit of a statistically isotropic distribution of points in S 2 as the 
number of points go to infinity. One can thus quantify anisotropy by calculating the departure of 
the mean observed probability distribution (Q bs(ci) ) from it, namely (Q bs{c>i) ) — &exp(cti) . 
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In practice, the expected <& e xp{on) for a statistically isotropic map is obtained sim- 
ply by scrambling a CMB map multiple times, and averaging over the resulting 
histograms. 

Lastly, to quantify anisotropy, we distill the histogram T(a^) into a single num- 
ber, by defining the indicator a = <r{9, (f>) as the variance of T(a^) (which has zero 
mean), namely 



-, N bin3 
1 y bins 



= 1 



Calculating a for the whole celestial sphere would yield a global measure of 
anisotropy. To obtain a directional indicator, we can instead calculate both MPASH 
and EPASH for spherical-shaped caps. The indicator a can then be viewed as a (non- 
local) measure of the anisotropy in the direction of the center of the cap. Thus, by 
construction, a(9, </>) measures the deviation from isotropy in a given direction, i.e. 
how the observed distribution of points deviates from a statistically isotropic one|B 

Since a is a discrete scalar function defined on S 2 , it can be expanded in spherical 
harmonics, and its power spectrum Di can be calculated, namely 

oo £ ^ 

a{9, <t>) = blm F * m ^> & and D * = 2£ + T ^ |5 ' m|2 ' (6) 

It then follows that if a large-angle asymmetry is present in the CMB temperature 
distribution, it should significantly affect the a— map on the corresponding angular 
scales (low-^ multipoles). 

In the next section, we shall generate the a— maps from LILC map with a Kp2 
sky cut, study its main fea ture s, and make a comparison with our previous results 
for the full-sky CMB mapsP^l 



3. Main Results and Conclusions 

Given that the large-scale angular correlations are nonlocal, a(8, <fi), calculated over 
a 30°-radius cap centered at (9, (/>), can be though of as a measure of the anisotropy 
in the direction (9, (f>). In our previous worlPZIthe strategy was to obtain a for a set 
of 12,288 caps of radius 30° co-centered with the same number of pixels generated 
by HEALPix with N S id e = 32, evenly covering the entire celestial sphere. The 
resulting directional map of anisotropy was the so-called a— map. We applied this 
new anisotropy indicator to three CMB WMAP maps: the LILCpSI and the TOH^ 
maps (which are two differently foreground cleaned full-sky maps resulting from the 
combination of the five frequency bands: K, Ka, Q, V, and W CMB maps measured 



b Incidentally, for a homogeneous and continuous distribution of points on S 2 the EPASH for a cap 
of aperture 9o < n/2 can also be calculated in a closed formp^l but in practice one approximates 
^ C exp ( a i) by tne MPASH ( & C e tp{ a i) ) °f tne statistically isotropic distribution of pixels. 
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by the WMAP satellite), and the co-added map, which is a weighted combination 
of the Q, V, and W WMAP maps. 

The resulting a— map were found to be anisotropic. Briefly, there is a prominent 
spot with very high a on the southeastern corner, with a well defined maximum 
at (b ~ 115°, I ~ 235°), which is close (by 16°) to the direction recently indicated 
in Ref. [11] It was further shown (by a standard spherical harmonics expansion) 
that the LILC a— map deviates from isotropy in a statistically significant way, with 
anomalously high (> 95% CL) dipole, quadrupole and octupole components (see 
Figs. 3 and 4 of Ref. |22|) . The higher components on the other hand fall within 
the expected values. This clearly indicates that the LILC map is not statistically 
isotropic. Finally, we noted that the quadrupole component has a very peculiar 
shape, being very symmetric around an axis slightly off the galactic North-South. 
Indeed, 82% of the total power in D2 comes from an axisymmetric component in 
the direction (b = 10°, I = 289°), somewhat close to the axes of symmetry of the 
temperature quadrupole and octupole found in Ref. S (about 24° from both) . 

As previously mentioned, however, there is no consensus as to whether the full 
sky cleaned maps available are indeed free of significant galactic contamination. The 
question then arises as to whether one should study the full sky maps or confine 
the analysis to regions where such contamination is small. 

In view of the lack of consensus on how to perform the data analysis, here we 
examine the robustness of our previous results by investigating the LILC map after 
the application of the Kp2 mask (hereafter the LILC-Kp2 map) , which discards the 
temperature fluctuations of 15.3% of the total number of pixels, mainly concentrated 
around the galactic plane. 

Note that some of the caps with centers close to but outside the Kp2 mask would 
still overlap with the mask itself. If the intersection region is too large the a value 
would be largely an artifact os the masking procedure. On the other hand, if we 
were to exclude all caps were any overlapping occurs, we would lose information on 
over half of the sky. To achieve a balance, we shall disregard caps that obey either 
of the following criteria: 

• have the cap center within the Kp2 mask, and 

• have over 15% of pixels within the Kp2 mask. 

With this critical value of 15% for the maximum number of overlapping pixels 
typically the value of a calculated for the same cap in both the LILC and LILC- 
Kp2 maps differ by less than 10%. 

As can be seen in Fig. [TJ there is a spot of very high a in the southeastern corner 
of the map, which coi nci des in both direction and magnitude with the one found in 
the full LILC a— map The fact that this spot lies well outside the region of signif- 
icant galactic contamination suggests it is not the result of galactic contamination. 
The possibility remains, however, that it is caused by some unaccounted foreground 
contamination. Such contamination however would unlikely to affect the different 
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frequency bands in exactly the same way. To verify whether this is the case, we 
calculated the a— map for the Q, V and W bands separately, along wi th the co- 
added map^3, which is considered the most reliable map for CMB studieJ^HSl (see 
Fig. [2]). The resulting maps are almost identical, supporting the result that fore- 
ground contamination may not account for the previously i 
anisotropy. 




0.020 0.30 



Fig. 1. The tr— map for the LILC-Kp2 map. This result was obtained by calculating a for spherical 
caps with aperture 8q = 30° , following the criteria outlined above. 



Q v 




Fig. 2. The cr— map for the WMAP's Q-, V- and W- bands, along with the combined co-added 
map, using the Kp2 mask. In all maps the high cr value is apparent. 

Several new smaller high a spots are also in evidence near the mask region, but 
they are probably just artifacts of the masking procedure. As discussed iiP^I the 
main features of the resulting cr— map are robust with respect to the number of 
spherical caps used to cover the celestial sphere or the cap aperture 9q. 
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To obtain more quantitative information about the observed anisotropy, we fol- 
lowed the procedure of our previous work and calculated the pow er spectrum of the 
LILC-Kp2 (j — map using the Anafast subroutine in Healpix-^ 1 . Since we are now 
dealing with an incomplete sphere, the spherical harmonics are no longer orthogo- 
nal, and the values obtained must be handled with care. The Di values are depicted 
in Fig. [3l along with the corresponding full sky values for comparison. It is clear 
that the dipole component of the cr— map is even larger for the LILC-Kp2 than for 
the full sky LILC. This is consistent with the presence of an axis of asymmetry, 
again confirming our earlier results. The direction of the a— map dipole changes, 
however, from (b = 141°, I = 240°) in the full sky map to (b = 150°, I = 209°), in 
the LILC-Kp2 map, a difference of 19°. The quadrupole and octupole components 
on the other hand are comparatively smaller in the LILC-Kp2 map, probably due 
to the fact that many of the high a structures other than the large spot in the 
southeastern quadrant are now excluded by the Kp2 mask. 



Fig. 3. The power spectrum of the LILC-Kp2 (black dots) and LILC (grey dots) cr— maps for 
1 = 1,..., 10. For comparison, the average cr— map power spectrum for a set of statistically isotropic 
full skies is shown (solid curve within the grey band), along with its 95% confidence limits (dashed 
curves) . 

Another feature missing from the LILC-Kp2 cr— map is the peculiar shape of the 
quadrupole component observed in the full sky map. In the latter case, 82% of the 
total power in D2 comes from an axisymmetric component in the direction (b = 
10°, I — 289°), which may be related to the observed alignment of the temperature 
quadrupole and octupole. 

We have shown that a large scale anisotropy, roughly resembling an axis of 
asymmetry, remains a feature of the CMB sky even if the regions where galactic 
contamination is large are not taken into account. Its essential dipolar nature in 
particular is a very robust feature, and is consistent with other axes reported in the 
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literature obtained using different methods. This strongly suggests that either the 
observable universe is intrinsically anisotropic, or that there are other, subtler forms 
of foreground contamination that have not yet been taken into account. Among the 
proposed explanations for the global preferred direction, it has been sugge sted t hat 
it could be due to a non-trivial topology of the spatial section of the universe^^l (f or 
more details on cosmic topology see the review articles Refs. l29l and, e.g., Refs. l30l 
Effl) . If topology is indeed the origin, the indicator T is promising in distinguishing 
between different topologies, as has been demonstrated by computer simulations in 
Ref.[531 These are very exciting possibilities, and are worthy of further investigation. 
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